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We analyze the possibility of encountering resonant transitions of high energy Majorana neu- 
trinos produced in Active Galactic Nuclei (AGN). We consider gravitational, electromagnetic and 
matter effects and show that the latter are ignorable. Resonant oscillations due to the gravitational 
interactions are shown to occur at energies in the PeV range for magnetic moments in the 10 _17 /is 
range. Coherent precession will dominate for larger magnetic moments. The alllowed regions for 
gravitational resonant transitions are obtained. 
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Majorana particles are natural representations of massive neutrinos since the most general mass term for a four 
component fermion field describes two Majorana particles with different masses. Majorana neutrinos also appear in 
many extensions of the minimal Standard Model; this is the case, for example, in SO(10) grand unified theories [1]. 

Neutrinos in general, and in particular Majorana neutrinos, can be used to probe the core of some of the most 
interesting cosmological objects. Due to their small cross sections these particles can stream out unaffected from 
even the most violent environments such as those present in Active Galactic Nuclei (AGN). The presence of several 
^>>0 ■ neutrino flavors and spin states modifies this picture: in their trek from their source to the detector the neutrinos 
can undergo flavor and/or spin transitions which can obscure some of the features of the source. Because of this, and 
due to the recent interest in neutrino astronomy (e.g. AMANDA, Nestor, Baikal, etc. [2]), it becomes important to 
understand the manner in which these flavor-spin transitions can occur, in the hope of disentangling these effects from 
the ones produced by the properties of the source. Without such an understanding it will be impossible to determine 
the properties of the AGN core using solely the neutrino flux received on Earth. 

In a previous publication [3] we considered the effects of the AGN environment on the neutrino flux under the 
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I ' assumption that all neutrinos were of the Dirac type. In this complementary publication we will consider the case of 
Oh, Majorana neutrinos and provide a deeper phenomenological study of this system, concentrating on the dependence 
of the effects on the magnitude of the neutrino magnetic moment and on the energy dependence of the predicted 
neutrino fluxes. 

The evolution of Majorana neutrinos in AGN is influenced both by its gravitational and electromagnetic interactions. 
The latter are due to the coupling with the magnetic field through a transition magnetic moment 1 . The combination 
of these effects leads to — ► z7 e or v r transitions. For simplicity we will deal with two neutrino species only, 
the extensions to three (or more) species is straightforward (though the analysis can become considerably more 
complicated [5]). 

The paper is organized as follows. We start with a brief description of neutrino production and gravitational 
oscillation in AGN environment in section 2. This is followed by a calculation of transition and survival probabilities 
of oscillating neutrinos and the resulting flux modifications (sections 3 and 4). In section 5 we give our conclusions. 
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II. NEUTRINO OSCILLATIONS IN AGN 



Active Galactic Nuclei (AGN) are the most luminous objects in the Universe, their luminosities ranging from 10 42 
to 10 48 ergs/sec. They are believed to be powered by a central supermassive black hole whose masses are of the order 
of 10 4 to 10 10 M Q . 



1 Majorana neutrinos are self conjugate particles which implies that the flavor diagonal magnetic moment form factor vanishes 
identically [4]. The absence of flavor-diagonal magnetic moments reduces the number of possible flavor and/or spin resonances 
that can be present compared to the Dirac case. 
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High energy neutrino production in an AGN environment can be described using the so-called spherical accretion 
model [6] from which we can estimate the matter density and also the magnetic field (both of which are needed to 
study the evolution of the neutrino system). Neutrino production in this model occurs via the ^ — ► ^ — > e ± 
decay chain [6,7] with the pions being produced through the collision of fast protons (accelerated through first-order 
diffusive Fermi mechanism at the shock [6]) and the photons of the dense ambient radiation field. These neutrinos 
are expected to dominate the neutrino sky at energies of 1 TeV and beyond [7]. 

Within this model the order of magnitude of the matter density p for typical cases can be estimated at 10 1 — 10 4 eV 4 . 
The magnetic field, for reasonable parameters, is of the order of 10 4 G [3,7,8]. 

In order to determine the effective interactions of the Majorana neutrinos in an AGN environment we start, follow- 
ing [3], from the Dirac equation in curved space including their weak and electromagnetic interactions 2 

[<7 a (^ + w„) - m+ / 75 + M<7 o6 F a6 ]V = 0, (1) 

where are the tetrads, j a the usual Gamma matrices, m is the mass matrix, / = J a -f a denotes the weak interaction 
current matrix, /j, is the neutrino magnetic moment matrix, F ab the electromagnetic field tensor, a ah — |[7 a ,7&]; and 
the spin connection equals 

^ = lha,-y b }e va ei !fl (2) 

where the semicolon denotes a covariant derivative. We used Greek indices (fi, v, . . .) to denote space-time directions, 
and Latin indices (a, b, . . .) to denote directions in a local Lorentzian frame. 

The method of extracting the effective neutrino Hamiltonian from (1) is studied in detail in [3]. We will therefore 
provide only a brief description of the procedure for completeness. 

The first step of the semiclassical approximation is to consider a classical geodesic x^(l) parameterized by an affine 
parameter Along this curve we construct three vector fields Vj^l), A = 1,2,3 such that x + i>a£, A satisfies the 
geodesic differential equation to first order in the t; A . We then use {l,£, A } as our coordinates (see, for example, [9]). 

Next we consider the classical action as a function of the coordinates which satisfies the relation p M = dS/ dx^ (p 
is the classical momentum), and define a spinor \ via the usual semiclassical relation 

V = e lS X (3) 

In our calculations it proves convenient to define a time-like vector p M corresponding to the component of momentum 
p orthogonal to the v^ A 

p " = p" - c A (A-yX; N AB = v\ V » B g^ (4) 

It can be shown that ca are constants [3]. We denote by R the length-scale of the metric, so that, for example, 
w M <~ 1/i?; and let P be the order of magnitude of the momentum of the neutrinos. We then make a double expansion 
of x, first in powers of £ and then in powers of 1/pR. We substitute these expressions into the Dirac equation and 
demand that each term vanish separately. 

It is then possible to reduce the resulting equations to a Schrodinger-like equation involving only \ which reads 

iX = H ce x 

ffeff = -hj + m) (*7 a e^+ H) + \™ 2 - —j+-t A l A + litO ah F ah , 
2 2 2 ^ 2 



(5) 



where, 



U = \NAB T B - <- v P^ V Al° 



2 The vector current is zero for Majorana neutrinos; they interact with matter only through the axial vector current unlike 
Dirac neutrinos. 
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(the over-bar indicates that the quantities are evaluated on the geodesic: x = x(l)). Eqn. (5) determines the evolution 
of the 8-component (for two flavors) spinor \- We are interested, however only in the 4-components representing the 
spinors with positive momenta [3] {i.e. those directed from the source to the observer). Projecting H c r into this 
4-dimcnsional subspace yields the effective Hamiltonian H c s for the states of interest 



• • 1 2 T ( 1 

H cS =ia+-m + p ■ J e ff I q _ 1 



p 



B* 
B 



where B denotes the magnetic field and the effective current J c g is defined by 

ja ja ja 

J cff — J W ~ J G 

in terms of the weak-interaction current and the "gravitational current" 

ja _ l f abcd- Xf (f . 2P f Pb \ 
J G - 7 e Afcd \ V b + — — J , 



where 



A/cd = {e f p, v - e fv ^)e%e v d , rf h = diag{\, -1, -1, -1). 



(6) 



(7) 



(8) 



(9) 



The term a in (6) is flavor and spin diagonal and can be eliminated by redefining the overall phase of \ with no 
observable consequences. 

The evolution of Majorana neutrinos through matter in the presence of strong gravitational fields incorporating 
magnetic effects is thus, 



(v e \ 








= H cS 


<z\ 













(10) 



where I is the affine parameter 3 and H c s is the 4x4 matrix containing the effects of the weak, electromagnetic and 
gravitational neutrino interactions, explicitly 



H e fi = 



( V ■ Jeff 



|Am 2 sin2i? 



V E/itB 



I Am 2 sin 2d 







p ■ Jeff + \ Am 2 cos 2tf -EntB* 

~Efl t B -p ■ Jeff 







Efi t B* 



\Am 2 sin2i? 



\ Am 2 sin 2d -p ■ J cB + \Am 2 cos 2-d ) 



(11) 



where ■& is the neutrino mixing angle, Am 2 — m\ — m\, E the energy of the particle, B = B\ +iB 2 where B\^ are the 
AGN magnetic field components perpendicular to the direction of motion, and \xt the transition magnetic moment. 
It is worth noting that, in contrast to the Dirac case, antincutrinos exhibit matter interactions. This Hamiltonian 
includes gravitational as well as electroweak effects; the latter have been studied previously (see, for example, [10]). 

In our calculations we will use the Kerr metric to allow for the possibility of rotation of the central AGN black hole 
(we also have assumed that the accreting matter generates a small perturbation of the gravitational field) . The metric 
for a Kerr black hole contains two parameters, r g , the horizon radius and a the total angular momentum of the black 
hole per unit mass. The geodesies in this gravitational field have three constants of the motion, commonly denoted by 



3 Note that the left hand side involves differentiation with respect to the affine parameter which has units of (mass) -2 . 
Therefore H e a has units of (mass) 2 which differs from the usual Hamiltonian units. In cases where the neutrino energy E is 
conserved the usual effective Hamiltonian is i7 e ff /E. 
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E, L and K. The first corresponds to the energy, the second to the angular momentum along the black- hole rotation 
axis; the third constant K has no direct interpretation, but is associated with the total angular momentum [11]. 

Using the AGN models mentioned above we can compare the weak interaction current Jw, to Jg, its gravitational 
counterpart. The orders of magnitude are, 



Jy 



w 



Gf P 



10 



-33 



in,, 



Pev) eV~ 



(12) 



where Gf is the Fermi coupling constant, m p the proton rest mass and p e y is the density in eV units which for 
typical cases is 10 1 — 10 4 [7]. Taking R ~ r g , the gravitational current part is found to dominate the weak current 
part for all relevant values of r g (10 14 to 10 20 eV~ 1 ). In the following we will drop Jw- 
Setting Jyy = 0, J c g can be written in terms of a dimcnsionlcss function /, 



\p- Jcft | = Er g 1 f(r/r g ,9,j,k,a/r g ). 
where we have chosen normalized parameters 



L 



K 



3 Er g {Er g f 



(13) 



(14) 



We have plotted \p ■ J e ff\ in figures 1 and 2 for some typical parameter values. 

Using (11) we can determine the AGN regions where resonant transitions occur. These resonances are governed 
by the 2x2 submatrices of (11) for each pair of states. The two possible resonances 4 are obtained by equating the 
diagonal terms for each submatrix and give rise to the following resonance conditions 



|Am 2 cos 2$ = p ■ J G 
i Am 2 cos 2$ = -p ■ J G 



(15) 



As can clearly be seen from the above two equations the resonant transitions do not occur simultaneously. We have 
considered electron and muon neutrinos, similar results hold for any other pair of flavors. 



log 1/1 




FIG. 1. Plot of In |/|, see (13), for j = 0.15, k = 0.4, a/r g = 0.4 and M = 10 M, 



4 The two remaining transitions have matrices with vanishing off-diagonal elements and will not exhibit resonances. 
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FIG. 2. Contour-plot of /, see (13), forj = 0.15, k = 0.4, a/r g = 0.4 for the case M = 1O 8 M . The contours a,b,c,d,e 
correspond respectively to / = ±1, ±10 _1 , ±10~ 2 , ±10~ 3 , ±10~ 4 , (positive values for < < tt/2, negative values for 
n/2<6< tt). 

Resonances occur provided / is comparable to ±Am 2 cos2dr g /E as can be seen from (15) and (16). As suggested 
by figures 1 and 2, we have verified that, just as in the Dirac case [3], for almost all values of j, k and a the neutrinos 
will undergo resonances at a radius r (the precise value of which changes with 9) for all relevant values of Am 2 (in 
the range 10~ 2 to lCT 10 eV 2 ) provided E is large enough (> 10 8 eV in this case) . This implies that practically all 
Majorana neutrinos will experience resonances provided their energy is large enough. As an example for E ~ 1 PeV 
and Am\ 2 ~ 10~ 6 eV 2 (Solar large angle solution), the resonance contour corresponds to (d) in Fig. 2. 

III. PROBABILITIES FOR ALLOWED TRANSITIONS 

In this section we evaluate the survival and transition probabilities of neutrino transitions and look into the region 
in parameter space where resonant neutrino transitions occur. The average probabilities for oscillating neutrinos 
(including non-adiabatic effects) produced in a region with mixing angle and detected in vacuum where the 
mixing angle is i9, are given in general by [12,13] 



2 - p LZ ) cos2t? G cos2tf 



(16) 



and 



where i)g is the gravitational mixing angle, 

tan 2d G 



P(v» -^Vc) = 7 :-[ 7 :-Plz) cos2i? g cos 2^ 



4E/j, t B 



Am 2 cos 2i? + 4p ■ J G 



(17) 



(18) 



prod. 



(evaluated at the production point) and \x t is the transition magnetic moment, B the magnetic field, Am 2 the usual 
mass difference parameter and Jq given in eqn.(13). Plz is the Landau Zener probability 



Plz = cxp <^ -27T 



a 



(19) 



where 



= E^tB\ 



a= dl [ - p - Jcs] < 



(20) 



The condition for adiabatic resonances to induce an appreciable transition probability is 2w 2 f3 2 > a which in terms 
of the magnetic moment implies 
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EB 



Am 2 



2tt 2 A 



1/2 



u ICS 

rmin 



(21) 



where A = \dl / din (p ■ J e s)\ is a measure of the scale of the gravitational field (divided by E), and we have assumed 
that the magnetic field remains constant over an interval of magnitude ~ A. In order to exhibit the energy dependence 
of /x™f n we estimate A ~ r g /E whence 



Mmin 



1 

B 



Am 2 



2ir 2 r g E 



1/2 



(22) 



Using equations (16) and (17) we can then calculate the survival and transition probabilities for the neutrinos. The 
results depend on the particular geodesic followed by the neutrinos, that is, it depends on the parameters E as well 
as j and k defined in (14). The results are also dependent on the characteristics of the gravitational field and its 
environment through the magnetic field, the horizon radius r g and the angular momentum parameter a. Finally there 
is also an important dependence on the neutrino parameters Am 2 and jji t . 

For the study of the probabilities there are two regions of interest (assuming (21) is satisfied). If p ■ Jq > (E (itB) 
the system experiences adiabatic resonances whenever the energy is such that da ~ 7r/2. li p ■ Jq < {E p t B) the 
system will exhibit coherent precession [14] if E /i t B > Am 2 or no appreciable transitions if E [i t B < Am 2 . 

For a study of the resonance scenario we will choose /it = 1CP 17 /ib which is amply allowed by the current ex- 
perimental and astrophysical constraints and still ensures neutrino adiabatic resonant transitions for all Am 2 values 
considered provided E is sufficiently large 5 . This value is smaller than the one usually considered in the literature 
due to our being concerned with energies in the PeV range and the fact that /U^L decreases as \j\fE. For larger 
values of p, t transitions are dominated by coherent precession (provided E[i t B > Am 2 ). 

Using expressions (16) and (17) we obtained the probabilities of gravitationally induced adiabatic transitions, the 
results are presented in Fig. 3. If E = En is the resonant energy (for a particular choice of parameters) then at that 
energy cos2$g ~ 0. H E < Er the transition probability approaches the vacuum value Piy^ — * v e ) = sin 2 2d. For 
E > En, and for the value of p, t chosen, -da — * t/2 and the transition probability approaches P(^ M — » D e ) = cos 2 
which is <~ 1; this is our region of interest. Equation (21) is satisfied in Fig. 3 when P(y^ — ► v e ) approaches 1 



and Ptyn 



approaches zero which corresponds to conditions for adiabatic resonant transitions. The energy 



En denotes the adiabatic threshold energy between no gravitational effects and adiabatic resonant conversion. This 
behavior is illustrated in Fig. 3. Fig. 4 shows the dependence of the probabilities on r and 9. 




18 20 
log E (eV) 




18 20 
log E (eV) 



FIG. 3. Plot of the persistence (left) and transition (right) probabilities as functions of E(eV) for /it = 10 17 /j-b, 
r/r g = 6, 8 = 7r/2, j = 0.15, k = 0.4, a/r g — 0.4 and M — 10 8 Mq. The curves a,b,c,d,e correspond to 
Am 2 = 10~ 10 , 10~ 8 , 10~ 6 , 10~ 4 , 10~ 2 eV 2 respectively. 



5 For E = IPeV and Am 2 = 10~ s eV 2 , A ~ r g /E, we have ~ 10~ 17 mb- 
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FIG. 4. Plot of the persistence (left) and transition (right) probabilities as functions of (r, 6) for [it = 10 17 fiB, j = 0.15, 
k = 0.4, a/r a = 0.4, E = IPeV, Am 2 = 10~ 8 eV 2 and M = 1O 8 M . 



The region in parameter space where adiabatic resonances occur is restricted by fi t > M^f n , \p ■ Jg\ > Efj, t B and 
by E > En (where En is the resonant energy corresponding to do = tt/4). Using (13) and (18) these conditions can 
be re-written as 



2 2(5 
" > "a, / > V, / > 5 

where 5 = r g Am 2 /4E and z/ = r g [itB. 

From (17) we obtain (for the conditions at hand Plz — 0) 



Pt(l-P t ) 



S + .f 




(23) 



(24) 



where P t = P{u^ — > v e ) denotes the transition probability and we have assumed d ~ 0. 
From (23) and (24) it follows that 



f + 6>f>v = (J 




(25) 



which required P t > (2 + V*2)/4 or P t < (2 - >/2)/4; since these expressions are invariant under the replacement 
P t -> 1 - P t we will assume P t > (2 + \/2)/4 ~ 0.85. 

For each fixed value of Pt the relations above define a region in the / — S — v space. The quantity / in (23) takes, 
for the situations under consideration, the values 10 — 1 depending on the production region, with the larger values 
associated with regions closer to the black hole horizon, (see Fig. 2). The allowed region in this space for two values 
of Pt are given in Fig. 5; the contour plot of v for various values of S and Pt are presented in Fig. 6 for a specific 
allowed value of / 



/ogv 




hgv 
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FIG. 5. Plot of v = r g it t B as functions of / and S = r g Am 2 /AE for transition probabilities P t = 0.87 (left) and P t = 0.99 
(right), We used M = 1O 8 M and B — 10 4 G and imposed (23). 



-P,=0.90 



-7.5 -7 -6.5 -6 



FIG. 6. Plot of v — r g fi t B as a function of 8 = r g Am 2 /AE for various values of the transition probability P t ; we chose 
M = 10 8 M Q , B = 10 4 G, / = 10~ 3 and imposed (23). 

As can be seen clearly from Fig. 6, for a reasonable value of / (1CP 3 ), adiabatic resonant transitions occur for fi t 
in the range ~ 1CT 18 - lCT 17 ^. 

The neutrinos will experience coherent precession if the conditions v > / and v > AS are satisfied; in this case 
the transition and persistence probabilities are ~ 1/2. For example, if / = 10~ 3 , coherent precession will occur fod 
fi t > 10 -1 Vb (for the above values of B and r g ) assuming that Am 2 < AE \i t B. For sufficiently large values of 
magnetic moment coherent precession will be the dominant mechanism. 



IV. FLUX MODIFICATION DUE TO GRAVITATIONAL OSCILLATIONS. 



According to the spherical accretion model [6] high energy neutrinos are produced by proton acceleration at the 
shock. Such acceleration is assumed to occur by the first order Fermi mechanism resulting in an E~ 2 spectrum 
extending up to £?£j OX . Estimates of expected neutrino fluxes from individual AGN are normalized, using this model, 
to their observed X-ray luminosities 6 ; the position of the maximum neutrino energy, however, is uncertain because 
it depends upon the parameters of the particular source. According to the model we have used [16], the neutrino 
flux at comparatively lower energies is strongly related to the X-ray flux; above 100 TeV the flux depends strongly on 
the turnover in the primary photon spectrum and one then needs to know the maximum proton energy E^^. The 
results of calculations of Szabo and Protheroe [16] give an approximate formula for the muon neutrino {y^ + v^) flux 
in terms of F x , the X-ray flux, and E^ ax 

F v ~ 0.25F X cxp(-20£/£P Q J£T 2 (26) 

where F x is a typical 2-10 Kev X-ray flux in erg cm _2 s _1 and E is the neutrino energy; the flux is plotted in (Fig. 7). 



6 There are varying results concerning the degree of proton confinement [15]. The various possibilities, however, alter only the 
low-energy neutrino spectrum and does not affect our results. 
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FIG. 7. Plot of the (f M + z/ M ) (solid) and {v e + v e ) (dashed) fluxes, as given by (26), as function of energy E when the X-ray 
flux equals F x = 10 _11 ergcm _2 sec _1 . E^ ax — 10 17 eV for AGN luminosity 10 45 ergs/sec, and : v e :: 2 : 1. 

Taking F x for a typical AGN to be 10 -11 ergcm _2 sec _1 (corresponding to for 3C273) and E^ iax = 10 17 eV for an 
AGN luminosity 7 L = 10 45 ergs/sec, the resulting neutrino flux is plotted in figure 5. 

As discussed in [17] and [18], the detectors AMANDA, Baikal and Nestor are sensitive to a wide range of neutrino 
parameters and will be able to test a variety of models of neutrino production in the AGN. For neutrino energies 
above 1 TeV, measuring the ultra high energy muon flux permits an estimation of the + flux, for neutrino 
energies above 3 PeV there is significant contribution to the muon rate due to v e interaction with electrons due to the 
W-resonance effect. Also, as shown in [18], the r rate from double bang events can be used to measure the v T + D T 
flux at energies 1 PeV and beyond. 

We noted earlier that matter effects are negligible in the AGN environment, but that gravity-induced resonances 
could cause a modification of the neutrino flux of any given flavor. This effect could cause an oscillation to r neutrinos 
8 generating a significant flux of t neutrinos to which planned experiments will be sensitive in the PeV range [18]. 
The observed r neutrinos would arise due to gravitationally induced resonance oscillations at the AGN itself. 

For the case of two flavors a and (3 the observed neutrino flux for species a, F a , can be expressed in terms of the 
initially produced neutrino fluxes F° and F2 as 



In order to calculate the fluxes in terms of neutrino oscillations we have made a number of simplifying but reasonable 
assumptions. 

• The initial (production) fluxes are assumed to be in the ratio : v e : v T :: 2 : 1 : 0. If charged and neutral 
pions are produced in equal proportions simple counting leads to equal fluxes of + and photons. The flux 
of v e + v e equals half the flux of + v^. Ref. [17] discusses the various models AGN neutrino fluxes. 

• There are equal numbers of neutrinos and antineutrinos. 

• As discussed above and in [3] matter effects are ignored. Resonances occur through the interaction of gravita- 
tional and electromagnetic interactions (proportional to the transition magnetic moment). We concentrate on 
the case where resonances are present corresponding to small magnetic moments. 

• It is assumed that there is no large matter effects in the path from AGN to the earth, so MSW effect does not 
cause oscillation for neutrinos in transit either. Vacuum oscillations do occur for exceedingly small Am 2 (for a 
distance of <~ 100 Mpc and energy 1 TeV, Am 2 < 10 _19 eV 2 ). Also the intergalactic magnetic field (~ 10~ 6 G) 
may cause spin transitions, but this is again a small effects in the parameter range we consider. 

Under these circumstances gravitationally induced spin flavor oscillations are the only important process which 
could cause neutrino transitions resulting in a modification of the observed fluxes. Figures 8 show the decrease of /i 



7 Assuming equipartition of magnetic energy and radiation energy at shock, [15] E^^is roughly proportional to \[L . where 
L denotes the AGN luminosity. 
8 There is negligible r neutrino production in the AGN environment according to all standard models . 



F a = [l- P(v, 




(27) 
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neutrino flux and the corresponding increase in the r neutrino flux produced by resonant oscillations of (i neutrinos to 
r neutrinos; the effect increases with energy and is very effective in the PeV range where r neutrinos can be detected 
relatively easily. 




FIG. 8. Flux modification of fi neutrino (left) and r neutrino fluxes due to gravitational oscillations of the — v T system 
in AGN environment. Am 2 = lCT 8 eV 2 , fj, t = W' 17 fi B , r/r g = 6, 6 = tt/2, j = 0.15, k = 0.4, a/r g = 0.4 and M = 10 8 M Q . 

The choice of [it — W~ 17 [1b and Am 2 = 10 _8 eV 2 for energy 1 PeV falls into the allowed region for gravitational 
resonant oscillations (Fig. 5) satisfying conditions (23). 



V. CONCLUSIONS 



We found that ultra-high energy Majorana neutrinos emanating from AGN are strongly affected by gravitational 
and electromagnetic effects. For typical values of gravitational current (Jg), which depends on the allowed black 
hole and geodesic parameters, gravitational resonant oscillations are found to occur for energies in which future 
neutrino telescopes [18] will be sensitive provided [it < 10~ 17 [1b- This therefore causes a significant increase of the 
corresponding r neutrino flux. Gravitational resonant oscillations cause flavor/spin oscillation of such neutrinos where 
the transition neutrino magnetic moments are small (Fig. 5). For larger magnetic moments, (fit > 10~ 17 /ib) coherent 
precession will dominate provided Am 2 < E fi t B 9 . The increase in the r neutrino flux is, of course, accompanied by 
a corresponding decrease in the fi neutrino flux. The two cases discussed above can in principle be differentiated by 
comparing the r and fi neutrino fluxes which would be equal for the case of coherent precession but markedly different 
for the case of resonant oscillations (Fig. 7 ). However, transition efficiency is greatest for the minimum value of the 
magnetic moment which corresponds to gravitational resonant oscillations. 

We have restricted our calculations to the two neutrino flavor case. A complete study should include at least 
three flavors and possibly four (to consider the possibility of sterile neutrinos). In the case of solar neutrinos the 
presence of more flavors can significantly alter the predicted fluxes [13]. However for the present situation where the 
experimental information on the AGN neutrino flux is quite limited, it is sufficient to determine the various effects 
and their strengths by using a a two flavor mixing description as we have done in our analysis. 
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9 The values of magneticmoment correspond to / = 10 3 . For neutrinos produced in the immediate vicinity of the horizon 
/ ~ 1 and resonances occur for [it < 10~ 13 /is while coherent precession dominates above this value. 
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